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We discuss 2D systems with Ising symmetry and competing interactions at different scales. In 
the framework of the Renormalization Group, we study the effect of relevant quartic interactions. 
In addition to the usual constant interaction term, we analyze the effect of quadrupole interactions 
in the self consistent Hartree approximation. We show that in the case of repulsive quadrupole 
interaction, there is a first order phase transition to a stripe phase in agreement with the well 
known Brazovskii result. However, in the case of attractive quadrupole interactions there is an 
isotropic-nematic second order transition with higher critical temperature. 
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The interest in phases with complex translational 
and/or orientational order is growing up in the con- 
densed matter community in the last years. Systems 
with very different microscopic origin present phase tran- 
sitions from disorder isotropic, homogeneous phases to 
anisotropic and/or inhomogeneous ones. The physical 
origin of this behavior is in many cases associated with 
the competition between short ranged attractive and 
long ranged repulsive interactions The attractive 

part tends to form ordered phases or condensates, or in 
the case of conserved order parameter, tends to produce 
phase separation. However, a long ranged repulsion frus- 
trates this tendency, favoring the emergence of complex 
phases that break translational and/or rotational symme- 
try. Examples of these systems go from highly correlated 
quantum systems like quantum Hall samples 0, 0] and 
High Tc superconductors 0| to classical systems like fer- 
romagnetic films 0, Q , diblock copolymers Q and liquid 
crystals 13, Hi, to cite but a few. All these systems have 
a special regime where physical properties are dominated 
by an extended region in momentum space with large de- 
generacy. For instance, in Fermi liquids this degeneracy 
is related with the existence of a Fermi surface at low 
temperatures. In classical systems, the competing inter- 
actions lead to a shift of the dominant wave vector in the 
structure factor to a non-zero value. In nearly isotropic 
systems, the low energy degrees of freedom are kinemat- 
ically constrained to a thin spherical shell of radius fcp 
determined by interactions. Examples of scalar and vec- 
tor order parameters behaving this way were studied by 
Brazovskii in a seminal work [8| . He showed that in sys- 
tems with a spectrum of fluctuations dominated by a 
non-zero wave vector, there is a first order phase transi- 
tion at a finite temperature form an isotropic to a stripe 
phase, induced by field fluctuations. His prediction was 
experimentally conflrmed in the microphase separation 
transition in diblock copolymers [§]. It is also observed 
in Monte Carlo simulations of ultrathin magnetic films 
with perpendicular anisotropy [l^ and in the Coulomb 



frustrated ferromagnet where the first order tran- 
sition is clearly seen. However, theoretical and numeri- 
cal work on ferromagnetic films 0, [l^ and the classic 
KTHNY theory of two dimensional melting [l^ , predict 
a more complex phase diagram. For instance, it could be 
possible to melt the stripes into a nematic phase, where 
the translational order is lost, but orientational order re- 
mains. 

With this motivation we review in this letter the ef- 
fective low energy theory for systems with competing 
interactions under the perspective of the Renormaliza- 
tion Group (RG)[l3]- The presence of a new scale ko 
with a large momentum space allowed for fluctuations, 
completely changes the properties of the Gaussian fixed 
point. We found that, upon expanding the angular mo- 
mentum content of the interaction, the d = 2 system is 
characterized by an infinite number of relevant quartic 
coupling constants. This result resembles the Fermi liq- 
uid theory where there is an infinite number of marginal 
Landau parameters controlling the fixed point [isj . In 
particular, we found that the first nontrivial coupling 
after the usual local (j)'^ theory, represents an interac- 
tion between local quadrupole moments. If this inter- 
action is repulsive, the Brazovskii analysis remains the 
correct one, predicting a first order transition to an inho- 
mogeneous state. However, if this interaction is attrac- 
tive, there is a new instability describing a second order 
isotropic-nematic phase transition with critical tempera- 
ture higher than the melting transition of the Brazovskii 
model. We have characterized this d — 2 nematic critical 
point in the self consistent Hartree approximation. We 
have computed the critical temperature Tc as well as the 
critical exponents /3 — 1/2, 7=1. The conditions for the 
existence of the nematic phase and its critical properties 
are the main results of this letter. In the rest of the paper 
we sketch the analysis leading to these results. 

In general, the low temperature physics of 2d models 
with competing interactions and Ising symmetry is well 
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described by a coarse-grained Hamiltonian of the type: 




where ro{T) ~ (T — Tc), k — \k\ and fco = |fco| is a con- 
stant given by the nature of the competing interactions. 
j^<Pk = Jq^ d6 J^^^l^ dk k and A ^ y/mro is a cut-off 
where the expansion of the free energy up to quadratic 
order in the momentum makes sense. The "mass" m 
measures the curvature of the dispersion relation around 
the minimum fco and the ". . ." in eq.([T]) means higher 
order terms in (fc — fco)- The structure factor has a max- 
imum at fc = fco with a correlation length ^ ~ l/y/mrQ. 
Therefore, near criticality (ro — > 0) the physics is domi- 
nated by an annulus in momentum space with momenta 
fc ~ fco and width 2A. This situation is quite similar with 
fermionic systems at low temperature, where the role of 
fco is the Fermi momentum, and the reduction of phase 
space to a spherical shell centered at the Fermi momen- 
tum is ruled by the Pauli exclusion principle. In our 
case of interest, the microscopic physics is very different, 
but the effects of kinematical constraints on momentum 
space are equivalent. 

We would like to determine what kind of interaction 
terms are relevant to study the low energy physics of a 
system given by the Hamiltonian of eq.([T]). The method 
is similar to the RG for Fermi liquids developed in r ef. 19l| 
and already applied to the Brazovskii model in ref. [20| to 
study the first-order transition to the stripe phase. The 
standard procedure is to identify a scale transformation 



that leaves the Gaussian theory invariant, and then study 
the relevance of interactions in the scaling limit very near 
the circle fc = fco- With this aim, as usual, we define 
a small radial wave vector q = fc — fco, we reduce the 
cutoff to A/s, with s > 1, and integrate over rapid modes 
A/s < \q\ < A, then we rescale the fields and wave vectors 
to reestablish the same scale and compare the couplings. 
First, note that the kinetic term in eq. |T]) is invariant 
under the rescalings 

A' = A/s, (2) 
q = qs, (3) 
= s--'/^cj,iq'/s). (4) 

Some important comments are in order. These trans- 
formations are independent of the dimension of the mo- 
mentum space. This result is very different for systems 
without competing interactions where fco = 0. In this 
case, the fields would scale as — s^^'^+^^/^(^((;'/s), 

where d is the spatial dimension of the system. This fact 
leads to a completely different analysis in the case of com- 
peting interactions. In fact, the usual concepts of upper 
and lower critical dimensions will change in our case, due 
essentially to the degeneracy of the lowest energy mani- 
fold. With this scaling in mind we immediately conclude 
that the term ro|0p in cq. ([1]) is relevant as it should 
be, as this term controls criticality. Let us now analyze a 
generic quartic interaction (we don't analyze cubic terms 
in this article since we are interested in systems with Ising 
symmetry (f> —(f)'- 
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Changing variables qi = ki — fcoi,with i — 1 . . .4, and 
rescaling the momenta and fields following eqs. ([2]), 
([3]) and ([4]), we obtain at "tree level" ^sj 94) = 

s^M(qi/s, 52/5, 93/5, (74/s)[2l|. We immediately see that 
the constant term is relevant. Therefore, at this level 
of approximation it is enough to keep the quartic term 
replacing = in the expression for u{qi). It is im- 
portant to note a difference with the case fco = 0; even 
though the coupling m is a constant, in the sense that 
it does not depend on g, it still depends on the angles 
Oi of each fcoi, and then u = u{9i^92, 0^,64^). However, 
these angles are strongly constrained by kinematics. Mo- 
mentum conservation allows us to eliminate one of them 
in terms of the other three. Furthermore, as the mo- 
menta are constrained to move in a very narrow circular 
region of radius fco and width A, we can fix only two of 



them and the other two are automatically slaved. Finally, 
since the system is rotational invariant, the couplings can 
only depend on the difference between these two angles, 
u{9i,92,0i + TT, 6I2 + vr) = u{9i - 6*2) = u{9). Thus, the 
quartic interaction is represented, not by few constants, 
but by a continuous function of an angle. There is still 
a last constraint, due to the fact that u(fci, fc2, fca, fc4) in 
Eq. ([5]) is invariant under any permutation of the four 
indexes, implying that u{9) — u{9 -I- tt). Therefore, we 
can expand u{9) in a Fourier series, obtaining an infinite 
set of coupling constants representing angular momen- 
tum "channels" of the interaction: 



u{9) = uq + U2 cos(20) + U4 cos(46') -f . . 



(6) 



The first term uq leads to the usual 0^ theory consid- 
ered by Brazovskii in his model for the isotropic-smectic 
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transition. In ref. [20| a detailed RG analysis showed the 
evolution of the running coupling constants uq and far 
from criticality, T << T^. This analysis is justified be- 
cause it is in this region that the fluctuation induced first 
order transition occurs. In this letter we are interested 
in another regime, T ^ Tc, where the whole set of cou- 
plings {uo, U2, U4, . . .} have to be considered, in principle, 
at the same level. A detailed analysis of the RG flow will 
be presented elsewhere 



where we have chosen Qij = a (nifij 



■^Sij^ and we have 
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For simplicity, let us consider the effect of the first two 
terms. A convenient way of representing the Hamiltonian 
in terms of uq and U2 in real space is 

Hint = J d?x [uo (b\x) + U2 tr Q^j , (7) 

where 



g,,(x) = 0(f) (^V,Vj - </)(x) (8) 

can be easily recognized as the quadrupole moment of 
the density 0^ 

Qij is a traceless symmetric tensor, being a natural 
local order parameter for a phase with orientational or- 
der. In fact, a phase with (0(f)) = 0, and Qij = 
J dx{Qij{x)) 7^ is a homogeneous phase, with orien- 
tational order with the nematic symmetry 9 9 + ir. 
In the following, we will show that under certain condi- 
tions, this is the ordered phase produced at the onset of 
the instability vq ^ 0. 

To study this phase transition, we make a self consis- 
tent Hartree approximation and analyze the Hamiltonian 
eq. ([7]) in the same lines of Brazovskii's work 0]. This 
approximation is exact in an 0{N) model in the limit 
N ^ 00. In our case, it will present corrections. 
In this analysis, we should also consider corrections of 
the type ^trQ'^ with 7 > 0. Although this term do not 
enter the critical properties of the system at this level of 
approximation, it will be important to stabilize the low 
temperature phase. 

As usual, we replace in eq. ^ (f)^ ^ (jP < (jP > and 
tr Cf- tr {(?!) (VjVj - ^V2(5y) 0} {Q^j{x)), where the 
mean values will be determined self consistently. With 
this procedure we obtain a quadratic Hamiltonian in the 
Hartree approximation given in momentum space by 
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Hartree 



J2 



'^-^^{k)[l3-'C-\k))^{-k), (9) 



where the static structure factor C{k) is now given by 

T 



C{k) 



and 



r -I- — (fc - fco)2 - ak"^ cos(26')(u2 + 7^^) 



r = ro + Uo 



(Pk 
(2^)2 



C(fc), 



, (10) 



(11) 



absorbed unimportant numerical factors in the definition 
of Uq. 9 is the angle subtended by k with the director 
n. Note that the new quartic term explicitly introduces 
an anisotropy in the structure factor of eq. PH]) . The 
coupling has exactly the form of the second term in eq. 
([6]), as it should be. This new term is also responsible for 
a shift in the value of the dominant wave vector. From 
the definition of the nematic order parameter, eq. ([8|l. 
we find that the amplitude of Qij is given by 



1 

" = -2 



<Pk 
(2^ 



fc2 cos(26i) C{k) 



(12) 



which, together with pil)) and (fTT|) completes a set of 
equations to be solved self-consistently. 

After the exact integration over the angles, we have 
integrated the radial variable k in the limit where the 
correlation length ^ ~ 1/ {^mr) is much larger than the 
typical wave-length of the system I/Zcq. It can also easily 
be checked that a = is always a solution of eq. (fTS]) 
and we expect that at high temperatures this is the only 
possible solution at finite r. Therefore, upon the k inte- 
grations at leading order in ^Jrrvr/kQ, the result can be 
expanded in powers of o? jr. It is convenient to write 
the equations in terms of the adimensional parameters 
r — > rT, fco — > koVmT and rg t — c(l — Tc/T), ob- 
taining at leading order in (a^ /r), 



, 1 

r = T + uokom—=z + U \ — 
\/r \ r 



and 



a^Uiir, T)+a2{r, T)—+0 




(13) 



0, (14) 



where ai{r,T) = U2{8 + U2k"§m^T'^ /2r^/'^) and a2{r,T) ^ 
8r7 -I- fc^m^U2(7 -I- k^U2 / r^) / ^/r . 

From this result, it is clear that if U2 > 0, correspond- 
ing to repulsive quadrupolar interaction, ai(r,T) > 
and a2{r,T) > for all r and T. Therefore, in this 
case, the only possible solution is a = 0. The model 
then reduces to that of Brazovskii, and a careful study 
of eq. (Uni) follows the same lines of ref. [§]. In this 
case nothing happens at T = T^. The system shows 
a first order phase transition to a stripe phase, with 
a melting temperature ~ Tc/{1 + rc) < Tc where 
Tc = r(T = 0) = {uokom)"^/^ . However, in the case of 
attractive quadrupole interactions, U2 < 0, eqs. ([T3ll and 
(|14p have non trivial solutions. For high temperatures 
T > Tc the only possible solution is a = as antic- 
ipated. This represents a high temperature disordered 
homogeneous and isotropic phase. On the other hand, 
if T < Tc, a nematic phase emerges continuously with 
a - |r - Tc|i/2. From eqs. ^ and ^ and the condi- 
tion ai{rc,Tc) = 0, we can read the critical temperature 
to be Tc = 4/(mfcQ) yju^ju^. 
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We have also computed the nematic susceptibihty, by 
couphng the system to a small external field conjugate 
to the nematic order parameter. Considering for sim- 
plicity the orientation of the field in the same direction 
than the order parameter, we find for the nematic sus- 
ceptibility Xn ^ ^/{T —Tc). This confirms the second 
order nature of this transition, with critical exponents 
(3 = 1/2 and 7 = 1. One has to bare in mind that these 
critical exponents will certainly be modified be fluctua- 
tions. The reason is that in our Hartree approximation 
the order parameter Qij depends on (j) fluctuations. In 
other words, it is a Hartree approximation for the order 
parameter 0, however a mean field in Qij. We expect 
that upon improving this approximation, the isotropic- 
nematic transition should be of the Kosterlitz-Thoules 
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typela^. 

Summarizing, we have analyzed a general model of a 
scalar field theory with competing interactions in two di- 
mensions. This model is dominated by a small circular 
shell in momentum space, profoundly modifying the crit- 
ical properties of the corresponding model without com- 
petition (fco = 0)- In particular, we have shown that in 
the framework of the Renormalization Group, the gaus- 
sian fixed point is affected by an infinite set of relevant 
coupling constants that codify the angular momentum 
content of a particular microscopic interaction. We have 
analyzed the simplest model with two coupling constants, 
corresponding to consider quadrupole moment interac- 
tions. We have found that the repulsive case reduces to 
the well known Brazovskii model, with a fluctuation in- 
duced first order phase transition to a non-homogeneous 
stripe state. However, for attractive quadrupole interac- 
tions the phase diagram changes considerably. By apply- 
ing a self-consistent Hartree approximation in the fields 
we found a phase transition from a high temperature 
disordered phase to a nematic phase at lower temper- 
atures. The critical temperature is higher than the melt- 
ing temperature found by Brazovskii, opening the pos- 
sibility of having a more complex phase diagram, with 
an intermediate nematic phase and a possible first order 
nematic-stripe phase transition. Whether U2 is positive 
or negative in a real system depends on the microscopic 
interactions. Our results also highlight how to measure 
nematic order from data for the structure factor and pre- 
dict the temperature window between and Tc opening 
the way to look for this phase experimentally, in sys- 
tems like anisotropic thin film ferromagnets and diblock 
copolymers, among many others. 
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